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(7j  =  1/y  ~  52  for  5  »  A2±).  Experimental  data 
confirm  this  behavior  (Fig.  3B).  This  dependence 
also  explains  why  quantum  jumps  were  not  ob¬ 
served  in  previous  experiments  with  NV  centers 
performed  at  low  magnetic  fields  [similar  mag¬ 
netic  field-enabled  decoupling  of  nuclear  spin 
was  proposed  recently  for  alkaline  earth  metal 
ions  (10,  22)].  The  dominance  of  flip-flop  pro¬ 
cesses  is  also  visible  in  the  quantum  state  tra¬ 
jectory  of  the  nuclear  spin  shown  in  Fig.  3C  (top). 
Here,  jumps  obey  the  selection  rule  A mj  =  +1 
imposed  by  the  flip-flop  term  HA.  From  analyz¬ 
ing  the  whole  quantum  state  trajectory,  a  matrix 
showing  the  transition  probabilities  can  be 
obtained  (Fig.  3C,  bottom). 

Single-shot  measurement  of  a  single  nuclear 
spin  places  diamond  among  leading  quantum 
computer  technologies.  The  high  readout  fidelity 
(92%)  demonstrated  in  this  work  is  already  close 
to  the  threshold  for  enabling  error  correction  (23), 
although  the  experiments  were  carried  out  in  a 
moderate-strength  magnetic  field.  Even  though 
the  optical  excitation  induces  complex  dynamics 
in  the  NV  center  (including  passage  into  singlet 
electronic  state),  the  nuclear  spin  relaxation  rates 
are  defined  solely  by  electron-nuclear  flip-flop 
processes  induced  by  hyperfine  interaction.  There¬ 
fore,  we  expect  improvement  of  7j  by  two  orders 
of  magnitude  (reaching  seconds  under  illumina¬ 
tion)  when  a  magnetic  field  of  5  T  is  used.  This 
will  potentially  allow  readout  fidelities  compara¬ 
ble  with  that  achieved  for  single  ions  in  traps 
(24).  The  present  technique  can  be  applied  to 
multiqubit  quantum  registers  (5,  6,  25),  enabling 


tests  of  nonclassical  correlations.  Finally,  single¬ 
shot  measurements  open  new  perspectives  for 
solid-state  sensing  technologies.  Spins  in  diamond 
are  considered  to  be  among  the  promising  candi¬ 
dates  for  nanoscale  magnetic  field  sensing  (26, 27). 
Currently  their  performance  is  limited  by  photon 
shot  noise  (26):  “Digital”  QND  will  provide  im¬ 
provement  over  conventional  photon  counting  in 
the  case  of  short  acquisition  time.  This  requires 
that  the  electron  spin  state  used  for  magnetic  field 
sensing  can  be  mapped  onto  the  nuclear  spin  with 
high  accuracy,  but  this  was  already  shown  to  be 
practical  in  NV  diamond  (5). 
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Strain-Induced  Pseudo-Magnetic 
Fields  Greater  Than  300  Tesla  in 
Graphene  Nanobubbles 

N.  I_evy,1,2* *t  S.  A.  Burke,1*^  K.  L.  Meaker,1  M.  Panlasigui,1  A.  Zettl,1'2  F.  Guinea,3 
A.  H.  Castro  Neto,4  M.  F.  Crommie1,2§ 

Recent  theoretical  proposals  suggest  that  strain  can  be  used  to  engineer  graphene  electronic  states 
through  the  creation  of  a  pseudo-magnetic  field.  This  effect  is  unique  to  graphene  because  of  its 
massless  Dirac  fermion-like  band  structure  and  particular  lattice  symmetry  (C3v).  Here,  we  present 
experimental  spectroscopic  measurements  by  scanning  tunneling  microscopy  of  highly  strained 
nanobubbles  that  form  when  graphene  is  grown  on  a  platinum  (111)  surface.  The  nanobubbles 
exhibit  Landau  levels  that  form  in  the  presence  of  strain-induced  pseudo-magnetic  fields  greater 
than  300  tesla.  This  demonstration  of  enormous  pseudo-magnetic  fields  opens  the  door  to  both 
the  study  of  charge  carriers  in  previously  inaccessible  high  magnetic  field  regimes  and  deliberate 
mechanical  control  over  electronic  structure  in  graphene  or  so-called  "strain  engineering." 


Graphene,  a  single  atomic  layer  of  carbon, 
displays  remarkable  electronic  and  me¬ 
chanical  properties  (1,  2).  Many  of  gra¬ 
phene’s  distinctive  properties  arise  from  a  linear 
band  dispersion  at  low  carrier  energies  (3)  that 
leads  to  Dirac-like  behavior  within  the  two- 
dimensional  (2D)  honeycomb  lattice — charge 
carriers  travel  as  if  their  effective  mass  is  zero 


(1).  An  intriguing  recent  prediction  is  that  a  dis¬ 
tortion  of  the  graphene  lattice  should  create  large, 
nearly  uniform  pseudo-magnetic  fields  and  give 
rise  to  a  pseudo-quantum  Hall  effect  (4).  Where¬ 
as  an  elastic  strain  can  be  expected  to  induce  a 
shift  in  the  Dirac  point  energy  from  local  changes 
in  electron  density,  it  is  also  predicted  to  induce 
an  effective  vector  potential  that  arises  from 


changes  in  the  electron-hopping  amplitude  be¬ 
tween  carbon  atoms  (5).  This  strain-induced  gauge 
field  can  give  rise  to  large  pseudo-magnetic 
fields  (Bs)  for  appropriately  selected  geometries 
of  the  applied  strain  (1,  6).  In  such  situations,  the 
charge  carriers  in  graphene  are  expected  to  cir¬ 
culate  as  if  under  the  influence  of  an  applied  out- 
of-plane  magnetic  field  (7-10).  It  has  recently 
been  proposed  that  a  modest  strain  field  with 
triangular  symmetry  will  give  approximately  uni¬ 
form,  quantizing  Bs  upward  of  tens  of  tesla  (4). 

Here,  we  report  the  measurement  of  Landau 
levels  (LLs)  arising  from  giant  strain-induced 
pseudo-magnetic  fields  in  highly  strained  graphene 
nanobubbles  grown  on  the  Pt(lll)  surface.  Lan- 
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Fig.  1.  STM  images  and  STS  spectra 
taken  at  7.5  K.  (A)  Graphene  mono- 
layer  patch  on  Pt(lll)  with  four 
nanobubbles  at  the  graphene-Pt  bor¬ 
der  and  one  in  the  patch  interior. 
Unreacted  ethylene  molecules  and  a 
small  hexagonal  graphene  patch  can 
be  seen  in  the  lower  right  (/tUnneiing  = 
50  pA,  l/sampie  =  350  mV,  3D  z-scale 
enhanced  4.6x).  (Inset)  High-resolution 
image  of  a  graphene  nanobubble  showing  distorted  honeycomb  lattice  resulting  from  strain  in  the 


bubble  (/tunneling  =  50  pA,  Vsamp[e  =  200  mV,  maxz  =  1.6  nm,  3D  z-scale  enhanced  2x).  (B)  STS 
spectra  of  bare  Pt(lll),  flat  graphene  on  Pt(lll)  (shifted  upward  by  3  x  10-11  ohm-1),  and  the 
center  of  a  graphene  bubble  (shifted  upward  by  9  x  10-11  ohm-1).  l/mod  =  20  mV. 
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Fig.  2.  (A)  Sequence  of  eight  dl/dV  spectra  ( T  ~  7.5  K,  l/mod  = 
20  mV)  taken  in  a  line  across  a  graphene  nanobubble  shown  in 
the  image  in  (B).  Red  lines  are  data  with  quartic  background 
subtracted;  black  dotted  lines  are  Lorentzian  peak  fits  (center  of 
peaks  indicated  by  dots,  with  blue  dots  indicating  n  =  0). 
Vertical  dash-dot  lines  follow  the  energy  progression  of  each 
peak  order.  (C)  Normalized  peak  energy  versus  sgn(n)^/\n\  for 
peaks  observed  on  five  different  nanobubbles  follow  expected 
scaling  behavior  from  Eq.  1  (dashed  line). 


dau  quantization  of  the  electronic  spectrum  was 
observed  by  scanning  tunneling  microscopy  (STM), 
which  revealed  pseudo-magnetic  fields  in  excess 
of  300  T.  Such  enormous  strain-induced  pseudo- 
magnetic  fields  may  allow  the  electronic  proper¬ 
ties  of  graphene  to  be  controlled  through  various 
schemes  for  applying  strain  (11),  as  well  as  the 
exploration  of  new  high-field  physical  regimes. 

Strained  graphene  nanobubbles  were  created 
by  in  situ  growth  of  sub-monolayer  graphene 
films  in  ultrahigh  vacuum  on  a  clean  Pt(  111)  sur¬ 
face  (12)  in  order  to  avoid  external  contamination 
and  trapped  gases.  The  epitaxial  graphene  was 
grown  by  exposure  of  Pt(l  1 1)  to  ethylene  followed 
by  annealing  (13,  14).  Graphene  grown  on  Pt  is 
expected  to  be  minimally  coupled  to  the  substrate, 
compared  to  graphene  grown  on  other  catalytic 
metals  (15,  16).  A  Dirac-like  band  structure  is 
preserved  for  graphene  on  Pt(l  1 1),  as  verified  by  a 
recent  photoemission  study  (17).  An  STM  to¬ 


pograph  of  the  graphene/Pt(lll)  surface  prepared 
in  this  manner  (Fig.  1A)  reveals  a  flat  graphene 
patch  (partially  surrounded  by  Pt)  that  encom¬ 
passes  five  graphene  nanobubbles.  Graphene  na¬ 
nobubbles  frequently  appear  near  the  edges  of  a 
graphene  patch,  but  are  also  sometimes  observed 
in  the  center  of  flat  patches  or  near  the  boundaries 
between  patches  and  are  presumably  pinned  near 
these  locations  (Fig.  1A).  These  nanobubbles  are 
likely  related  to  the  larger-scale  “wrinkle”  struc¬ 
tures  observed  by  low  energy  electron  microscopy 
that  form  upon  cooling  as  a  result  of  the  differing 
thermal  expansion  coefficients  of  graphene  and  the 
platinum  surface  (1 7). 

Individual  nanobubbles  often  have  a  triangu¬ 
lar  shape  (Fig.  1A,  inset),  reflecting  the  lattice 
symmetry  of  the  graphene  and  the  underlying  Pt 
surface,  and  are  typically  4  to  10  nm  across  and 
0.3  to  2.0  nm  tall.  Atomic-resolution  imaging  of 
the  nanobubbles  confirms  the  honeycomb  struc¬ 


ture  of  graphene  here  (Fig.  1A,  inset),  although 
the  lattice  is  distorted  because  of  the  large  strain 
occurring  in  these  structures.  The  expected  strain- 
induced  pseudo-magnetic  field  in  a  graphene  na¬ 
nobubble  can  be  estimated  by  using  the  relation 
O  =  (p/z2//a)O0  for  the  flux  per  ripple  in  a  dis¬ 
torted  graphene  sheet  (6),  where  h  is  the  height,  / 
is  the  width,  a  is  on  the  order  of  the  C-C  bond 
length,  and  O0  is  the  quantum  of  flux.  The  param¬ 
eter  p  =  dlog(^)  /  Slog  (a)  relates  the  change  in  the 
hopping  amplitude  between  nearest  neighbor 
carbon  atoms  (t)  to  bond  length  and  has  a  typical 
magnitude  of  2  <  |5  <  3  for  graphene.  For  a 
nanobubble  of  /  =  4  nm  and  h  =  0.5  nm,  this  yields 
a  Bs  of  order  100  T.  The  large  curvature  and 
correspondingly  high  strain  incorporated  into 
the  triangular  nanoscale  bubbles  observed  here 
make  them  ideal  candidates  for  the  observation 
of  pseudo-LL  because  of  large  strain-induced 
pseudo-magnetic  fields. 

The  local  electronic  structure  of  strained  graphene 
nanobubbles  and  surrounding  graphene  films 
was  characterized  by  scanning  tunneling  spec¬ 
troscopy  (STS)  performed  at  -  7.5  K  by  using 
standard  lock-in  techniques  to  obtain  differen¬ 
tial  conductance  (dl/dV).  The  measurement  of 
dl/dV  reflects  features  in  the  local  density  of 
states  (LDOS)  of  the  surface  at  the  position  of  the 
STM  tip  (18).  STS  measurement  of  the  bare  Pt 
surface  was  used  to  calibrate  the  LDOS  of  the  tip 
upon  approach  and  between  sequences  of 
spectra.  STS  spectra  measured  over  the  bare  Pt 
regions  (Fig.  IB)  are  relatively  featureless  and 
show  the  expected  Pt(lll)  surface  state  (19). 
Spectra  recorded  over  the  flat  graphene  patches 
show  a  subtly  modified  structure  compared  with 
the  clean  Pt(  111)  surface,  and  no  clear  signatures 
of  the  graphene  Dirac  point  were  observed  in 
these  regions  (Fig.  IB).  Spectra  measured  at  the 
boundary  between  the  flat  graphene  and  the 
nanobubbles  (fig.  S2)  exhibit  features  consistent 
with  a  Dirac  point  located  -300  mV  above  the 
Fermi  energy,  as  recently  observed  by  photoemis¬ 
sion  (1 7),  as  well  as  a  gaplike  feature  with  a  full 
width  at  half  maximum  (FWHM)  of  127  ±  9  mV 
centered  at  the  Fermi  energy  (Fsampie  =  0)  recently 
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Fig.  3.  (A)  Experimental  topographic  line  scan  and  experimentally  determined 
Bs  profile  over  the  tip  trajectory  shown  by  black  line  in  (B).  (B)  STM  topography 
of  graphene  nanobubble.  (C)  Topography  of  theoretically  simulated  graphene 
nanobubble  with  calculated  Bs  color  map.  (D)  Simulated  topographic  line  scan 
and  Bs  profiles  extracted  from  line  shown  in  (C). 


associated  with  phonon-assisted  inelastic  tunnel¬ 
ing  (20). 

STS  measurements  made  directly  over  the 
nanobubble  regions  exhibit  a  succession  of  rela¬ 
tively  strong  peaks,  spaced  by  more  than  100 
meV,  which  do  not  appear  in  the  spectra  on  other 
regions  of  the  sample  (Fig.  IB).  The  peaks  are 
typically  weaker  at  negative  bias,  which  may  be 
attributed  to  the  expected  shorter  vertical  exten¬ 
sion  of  wavefunctions  at  lower  energies.  These 
peaks  in  the  LDOS  of  the  graphene  nanobubble 
are  unlike  features  seen  previously  in  STS  per¬ 
formed  on  graphene  in  the  absence  of  a  magnetic 
field  on  SiC  (21, 22)  and  Si02  (20, 23)  substrates 
and  often  overwhelm  the  usual  graphene  fea¬ 
tures.  Figure  2  A  shows  a  series  of  spectra  taken  at 
different  positions  across  a  single  nanobubble 
(topography  shown  in  Fig.  2B).  These  spectra 
display  the  typical  nanobubble  peak  structure  as 
well  as  the  inelastic  feature  at  the  Fermi  energy 
described  above  [a  broader  bias  range  is  shown  in 
(fig.  S3)].  Similar  peak  structure  was  observed  on 
10  different  nanobubbles  with  four  different  STM 
tips,  with  some  variation  in  peak  spacing  and 
amplitude  that  we  presume  was  caused  by  varia¬ 
tions  in  strain-induced  electronic  structure  arising 
from  different  nanobubble  geometries. 

These  peaks  observed  in  the  nanobubble 
LDOS  can  be  attributed  to  LLs  originating  from 
a  strain-induced  pseudo-magnetic  field.  Other 
possible  origins  that  we  mle  out  as  unlikely  in¬ 
clude  confinement  effects  and  defect  creation. 
Confinement  of  Dirac  fermions  in  graphene  is 
difficult  because  of  Klein  tunnelling  and  the  sup¬ 
pression  of  backscattering  (1).  If  confinement 
somehow  occurred  at  the  nanobubble  edges,  the 
peaks  observed  would  be  expected  to  follow  a 
different  progression,  most  notably  missing  the 
n  =  0  peak  observed  in  spectra  taken  at  the 
center  of  most  nanobubbles  (24).  Confinement 
would  also  result  in  strong  nodal  patterns  in 
LDOS  (24),  which  were  not  observed.  The  pro¬ 
duction  of  defect  states  by  strain  is  also  unlikely 
because  of  the  high  energetic  barrier  toward  de¬ 


fect  creation.  It  is  possible  that  some  nanobub¬ 
bles  were  formed  at  defect  sites  in  the  as-grown 
graphene,  but  the  nanobubbles  that  we  mea¬ 
sured  did  not  show  signatures  of  defect  physics 
in  spectroscopy,  such  as  tip  height  dependence 
of  peak  positions  associated  with  the  charging  of 
defect  states. 

The  most  likely  explanation  for  the  nano¬ 
bubble  peaks  is  that  they  arise  from  a  large, 
relatively  uniform  strain-induced  pseudo-magnetic 
field.  This  pseudo-magnetic  field  is  expected  to 
mimic  the  influence  of  a  real  magnetic  field  ap¬ 
plied  perpendicular  to  the  graphene  sheet  and 
give  rise  to  LLs  (4).  These  appear  as  a  series  of 
peaks  in  STS  as  they  do  for  the  case  of  a  real 
magnetic  field  (25,  26).  Specifically,  the  2D 
massless  Dirac  nature  of  charge  carriers  near  the 
Fermi  energy  causes  the  progression  of  peaks  in 
the  LDOS  to  follow  the  expression  (7,  8): 

En  =  sgn(n)hacy/\n\  +  E^c, 

C0c  =  yj2  eh^Bs  (1) 

where  En  is  the  position  in  energy  of  the  nth  LL 
with  respect  to  the  Fermi  level,  coc  is  the  cyclotron 
resonance  frequency  arising  from  Bs,  and  vF  is  the 
Fermi  velocity.  Both  positive  and  negative  n  will 
appear  symmetrically  about  the  Dirac  point, 
corresponding  to  electron  and  hole  states,  respec¬ 
tively,  as  well  as  mn  =  0  state  coincident  with  the 
Dirac  point  (1,  7,  8). 

The  energy  progression  expected  for  LLs  in 
graphene  can  be  compared  to  STS  data  taken  on 
a  nanobubble  by  fitting  the  spectra  with  a  se¬ 
quence  of  Lorentzian  peaks  following  Eq.  1  with 
a  simple  polynomial  background  (Fig.  2 A).  The 
observed  peak  structure  follows  the  expected 
progression  well,  and  a  value  of  Bs  for  position 
“1”  in  Fig.  2A  was  determined  to  be  350  ±  40  T 
(energies  of  the  n  =  0,  1,  and  2  states  from  the  fit 
are  shown  below  the  spectrum  for  this  position). 
Additional  peaks  that  largely  follow  the  expected 
sgn(w)  y/\n  \  progression  were  also  observed  over 


a  wider  bias  range.  Small  deviations  from  this 
progression  occur  for  higher  energy  peaks,  be¬ 
cause  the  graphene  dispersion  is  not  strictly  lin¬ 
ear  in  this  range.  A  plot  of  normalized  energy 
(En  - E Dh-ac) /  versus  sgn (n)  y/\n\,  shown  in 
Fig.  2C,  compiled  from  spectra  on  five  different 
nanobubbles,  demonstrates  the  expected  scaling 
behavior  for  LLs  in  graphene. 

The  offset  and  spacing  of  the  peak  progres¬ 
sion  changes  over  different  positions  on  individ¬ 
ual  nanobubbles,  indicating  a  spatial  variation 
of  Bs  and  £Dirac.  The  Bs  profile  across  a  strained 
nanobubble  was  extracted  from  the  spectral  peak 
spacing  at  different  positions  across  a  nanobub¬ 
ble,  shown  in  Fig.  3B.  As  seen  in  Fig.  3 A,  the  Bs 
profile  is  reasonably  flat  across  the  center  of  the 
bubble,  indicating  a  relatively  uniform  pseudo- 
magnetic  field  of  300  to  400  T  for  this  particular 
geometry.  EDkac,  coincident  with  the  n  =  0  LL, 
ranges  from  0.2  to  0.3  eV  across  this  region  of 
the  nanobubble.  This  variation  in  ED^rac  indicates 
scalar  potential  variations  across  the  nanobubble, 
as  expected  for  an  elastic  deformation  that  does 
not  constitute  a  pure  shear  strain  (1,  5). 

To  compare  the  experimental  spatial  depen¬ 
dence  of  the  strain-induced  pseudo-magnetic 
field  to  theoretical  predictions,  we  simulated  a 
triangular  nanobubble  with  similar  geometry  to 
that  shown  in  Figs.  2A  and  3B  by  using  con¬ 
tinuum  elasticity  theory  (27)  [see  supporting  in¬ 
formation  for  details  (12)].  In  this  calculation,  the 
edges  of  a  triangular  graphene  patch  were  brought 
in  toward  the  center  to  simulate  the  strain  arising 
from  the  different  coefficients  of  thermal  expan¬ 
sion  of  graphene  and  the  underlying  Pt(  111)  sub¬ 
strate  during  experimental  sample  preparation. 
The  predicted  pseudo-magnetic  field  arising 
from  this  strain  field  was  also  calculated  follow¬ 
ing  (4).  A  3D  plot  of  the  simulated  nanobubble 
shape  with  the  pseudo-magnetic  field  strength 
shown  as  a  color  map  is  displayed  in  Fig.  3  C  (the 
corresponding  experimental  nanobubble  is  shown 
in  Fig.  3B).  The  simulated  and  experimental  topo¬ 
graphic  profiles  agree  well,  as  seen  by  the  corn- 
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parison  between  Fig.  3,  A  and  D.  The  simulated  Bs 
profile  (Fig.  3D)  also  agrees  reasonably  well  with 
the  experimentally  extracted  Bs  (Fig.  3A)  given 
that  there  are  some  differences  in  shape  and  un¬ 
certainty  in  the  exact  boundary  conditions.  The 
main  discrepancy  is  that  the  experimental  Bs  pro¬ 
file  appears  much  more  uniform  than  the  sim¬ 
ulated  one,  which  may  be  attributed  to  spatial 
averaging  of  the  LDOS  over  the  magnetic  length 
scale  of  ~1.5  nm.  The  simulated  and  experimen¬ 
tal  nanobubbles  both  exhibit  a  relatively  uniform 
strain-induced  pseudo-magnetic  field  of  200  to 
400  T  across  the  central  region  with  increasing 
field  at  the  edges. 

The  exceptional  flexibility  and  strength  of  gra¬ 
phene  membranes  (28, 29)  coupled  with  the  large 
strain-induced  fields  observed  suggest  that  strain 
engineering  of  nanoscale  energy  levels  (11,  30) 
may  be  a  viable  means  of  controlling  the  elec¬ 
tronic  structure  of  graphene,  even  at  room  temper¬ 
ature.  The  experimental  demonstration  of  these 
enormous  pseudo-magnetic  fields  also  provides 
a  new  basis  for  the  study  of  extreme  high  magnet¬ 
ic  field  regimes  in  a  condensed-matter  environment. 
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Probing  the  Superfluid-to-Mott 
Insulator  Transition  at  the 
Single-Atom  Level 

W.  S.  Bakr/  A.  Peng/  M.  E.  Tai/  R.  Ma/  ].  Simon/  ].  I.  Gillen/  S.  Foiling/'2 
L.  Pollet/  M.  Greiner1* * 

Quantum  gases  in  optical  lattices  offer  an  opportunity  to  experimentally  realize  and  explore 
condensed  matter  models  in  a  clean,  tunable  system.  We  used  single  atom-single  lattice  site 
imaging  to  investigate  the  Bose-Hubbard  model  on  a  microscopic  level.  Our  technique  enables 
space-  and  time-resolved  characterization  of  the  number  statistics  across  the  superfluid-Mott 
insulator  quantum  phase  transition.  Site-resolved  probing  of  fluctuations  provides  us  with  a 
sensitive  local  thermometer,  allows  us  to  identify  microscopic  heterostructures  of  low-entropy 
Mott  domains,  and  enables  us  to  measure  local  quantum  dynamics,  revealing  surprisingly  fast 
transition  time  scales.  Our  results  may  serve  as  a  benchmark  for  theoretical  studies  of  quantum 
dynamics,  and  may  guide  the  engineering  of  low-entropy  phases  in  a  lattice. 


Microscopic  measurements  can  reveal 
properties  of  complex  systems  that  are 
not  accessible  through  statistical  en¬ 
semble  measurements.  For  example,  scanning 
tunneling  microscopy  has  allowed  physicists  to 
identity  the  importance  of  nanoscale  spatial  in- 
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homogeneities  in  high-temperature  superconduc¬ 
tivity  (1),  and  single-molecule  microscopy  (2)  has 
enabled  studies  of  local  dynamics  in  chemical 
reactions,  revealing  the  importance  of  multiple 
reaction  pathways  (3).  Whereas  previous  ultra¬ 
cold  quantum  gas  experiments  focused  primar¬ 
ily  on  statistical  ensemble  measurements,  the 
recently  introduced  single  atom-single  lattice  site 
imaging  technique  in  a  quantum  gas  microscope 
(4)  opens  the  door  for  probing  and  controlling 
quantum  gases  on  a  microscopic  level.  Here,  we 
present  a  microscopic  study  of  an  atom-lattice 


system  that  realizes  the  bosonic  Hubbard  model 
and  exhibits  a  quantum  phase  transition  from  a 
superfluid  to  a  Mott  insulator  (5-7). 

In  the  weakly  interacting  superfluid  regime, 
the  many-body  wave  function  factorizes  into  a 
product  of  states  with  well-defined  phase  on  each 
lattice  site,  known  as  coherent  states,  with  Pois- 
sonian  number  fluctuations.  As  the  strength  of 
the  interaction  increases,  the  number  distribution 
is  narrowed,  resulting  in  a  fixed  atom  number 
state  on  each  site  deep  in  the  Mott  insulator  re¬ 
gime.  We  study  this  change  in  the  number  sta¬ 
tistics  across  the  transition;  these  microscopic 
studies  are  complementary  to  previous  experi¬ 
ments  that  have  focused  on  measuring  ensemble 
properties  such  as  long-range  phase  coherence, 
excitation  spectra,  or  compressibility  (7-9).  Local 
properties  such  as  on-site  number  statistics  (10) 
had  been  accessible  only  indirectly  (8,  11,  12) 
and  were  averaged  over  several  shells  of  super¬ 
fluid  and  Mott  insulating  domains  in  the  inho¬ 
mogeneous  system,  complicating  quantitative 
interpretation.  More  recently,  the  shell  structure 
was  imaged  through  tomographic  (13),  spec¬ 
troscopic  (14),  and  in  situ  imaging  techniques, 
coarse-grained  over  several  lattice  sites  (15). 

We  started  with  a  two-dimensional  (2D)  87Rb 
Bose-Einstein  condensate  of  a  few  thousand  atoms 
confined  in  a  single  well  of  a  standing  wave,  with 
a  harmonic  oscillator  length  of  130  nm  (16).  The 
condensate  resided  9  pm  from  an  in-vacuum  lens 
that  was  part  of  an  imaging  system  with  a  resolu¬ 
tion  of  -600  nm.  This  high-resolution  system 
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